Field propagation in the Schwarzschild-de Sitter black hole 
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We present an exhaustive analysis of scalar, electromagnetic and gravitational perturbations in 
the background of a Schwarzchild-de Sitter spacetime. The field propagation is considered by means 
of a semi-analytical (WKB) approach and two numerical schemes: the characteristic and general 
initial value integrations. The results are compared near the extreme regime, and a unifying picture 
is established for the dynamics of different spin fields. Although some of the results just confirm 
usual expectations, a few surprises turn out to appear, as the dependence on the non-characteristic 
initial conditions of the non-vanishing asymptotic value for I = mode scalar fields. 



PACS numbers: 04.30.Nk,04.70.Bv 



I. INTRODUCTION 



Wave propagation around non trivial solutions of Ein- 
stein equations, black hole in particular, is an active field 
of research (see 0. Q and references therein) . The per- 
spective of gravitational waves detection in a near future 
and the great development of numerical general relativity 
have increased even more the activity on this field. Grav- 
itational waves should be especially strong when emitted 
by black holes. The study of the propagation of per- 
turbations around them is, hence, essential to provide 
templates for the gravitational waves identification. In 
the other hand, recent astrophysical observations indi- 
cate that the universe is undergoing an accelerated ex- 
pansion phase, suggesting the existence of a small posi- 
tive cosmological constant and that de Sitter (dS) geome- 
try provides a good description of very large scales of the 
universe 0. We notice also that string theory has re- 
cently motivated many works on asymptotically anti-de 
Sitter spacetimes (see, for instance, H, IE ) ■ 

In this work, we perform an exhaustive investigation of 
scalar, electromagnetic and gravitational perturbations 
in the background of a Schwarzchild-de Sitter spacetime. 
We scan the full range of the cosmological constant, from 
the asymptotic flat case (A = 0) up to the critical value 
of A which characterizes the Nariai solution Two 
different numerical methods and a higher order WKB 
analysis are used. 

We remind that for any perturbation in the space- 
times we consider, after the initial transient phase there 
are two main contributions to the resulting asymptotic 
wave initially the so called quasinormal modes, 

which are suppressed at later time by the tails. The first 



can be understood as candidates to normal modes which, 
however, decay (their energy eigenvalues becomes com- 
plex), as in the ingenious mechanism first described by G. 
Gamow in the context of nuclear physics ^2 ■ After the 
initial transient phase, the properties of resulting wave 
are more related to background spacetime rather than to 
the source itself. 

It is well known that for asymptotic flat backgrounds 
the tails decay according to a power-law, whereas in a 
space with a positive cosmological constant the decay is 
exponential. Curiously, I = modes for scalar fields in 
dS spacetimes, contrasting to the asymptotic flat cases, 
approach exponentially a non vanishing asymptotic value 
Ha Hi- We detected, by using a non-characteristic nu- 
merical integration scheme, a dependence of this asymp- 
totic value on the initial velocities. In particular, it van- 
ishes for static initial conditions. 
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The semi-analytical analyses of this work were per- 
formed by using the higher order WKB method proposed 
by Schutz and Will [lj|, and improved by Iyer and Will 
(la. Il7| . It provides a very accurate and systematic way 
to study black hole quasinormal modes. We apply it 
to the study of various perturbation fields in the non- 
asymptotically flat dS geometry. Quasinormal modes are 
also calculated according to this approximation, and the 
results are compared to the numerical ones whenever ap- 
propriate, providing a quite complete picture of the ques- 
tion of quasinormal perturbations for dS black holes. 

Two very recent works overlap our analysis presented 
here. A similar WKB analysis [3 wa s done very re- 
cently by Zhidenko [l9|. and his results coincide with 
ours. Yoshida and Futamase [2(j used a continued frac- 
tion numerical code to calculate quasinormal mode fre- 
quencies, with special emphasis to high order modes. Our 
results are also compatible. Finally, we notice that solu- 
tions of the wave equation in a non trivial background has 
also been used to infer intrinsic properties of the space- 
time E3. 
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II. METRIC, FIELDS AND EFFECTIVE 
POTENTIALS 



The metric describing a spherically black hole in the 
presence of a cosmological constant is well known in 
the literature. Written in spherical coordinates, the 
Schwarzschild-de Sitter metric is given by 

ds 2 = -h{r)dt 2 + h(r)- l dr 2 +r 2 (dO 2 + sin 2 9 d(f> 2 ) , (1) 



where the function h(r) is 



h(r) = 1 - 



2m 



Ar 2 
~3~ 



(2) 



The integration constant m is the black hole mass, and 
if the cosmological constant is positive, the spacetime is 
asymptotically de Sitter. In this case, A is usually written 
in terms of a "cosmological radius" a as A = 3/a 2 . 

Assuming m > and < 9m 2 A < 1, the function 
h(r) has two positive zeros r + and r c and a negative 
zero r„ = — (r+ +r c ). This is the Schwarzschild-de Sitter 
(SdS) geometry, in which we are interested. The horizons 
r + and r c are denoted event and cosmological horizons 
respectively. In this case, the constants m and a are 
related with the roots by 



2 2,2, 

a = r + + r c + r + r c 



2ma 



r+r c (r + + r c ) 



(3) 



(4) 



If 9m 2 A = 1, the zeros r + and r c degenerate in a double 
root. This is the extreme SdS black hole. If 9m 2 A > 1, 
there is no real positive zeros, and the metric (JIJ docs 
not describe a black hole. 

We will consider scalar, electromagnetic and gravita- 
tional perturbations in the submanifold given by 



T+ = {(t,r,6,, 



< r <r c } 



(5) 



In this region, we define the "tortoise" radial coordinate 
by 

x(r) = — — ln(r c - r) - — — ln(r - r+) + —— hx(r - r n ) , 

Zit\i£ Zitxifi 

(6) 



with 



dhir) 



dr 



(7) 



The constants n + and the surface gravity associ- 

ated with the event and cosmological horizons. 

For a scalar field $ obeying the massless Klein-Gordon 
equation 



□$ = , 



(8) 



the usual separation of variables in terms of a radial field 
and a spherical harmonic Y^ m (#, tp), 



(9) 



leads to Schrodinger-type equations in the tortoise coor- 
dinate for each value of £, 



d 2 ipf c d 2 ip 



2,/,sc 



+ 



dt 2 dx 2 
where the effective potential V sc is given by 



V sc (r) = h(r) 



£(£ + 1) 



2m 



(10) 



(11) 



In the Schwarzschild-de Sitter geometry, in contrast 
to the the case of a electrically charged black hole, it is 
possible to have pure electromagnetic and gravitational 
perturbations. For the first, the potential of the corre- 
sponding Schrodinger-type equation is given |22| 



V el (r) = h(r) 



£(£ + 1) 



(12) 



with I > 1. The gravitational perturbation theory for 
the exterior Schwarzschild-de Sitter geometry has been 
developed by 0, • The potential for the axial and polar 
modes are, respectively, 



V ax (r) = h(r) 



£{£+!) 6m 



(13) 



V*°(r) = 



2h(r) 
r 3 (3m + cr)' 



[9m 3 + 3c 2 mr 2 + c 2 (l + c)r 3 



+3m 2 (3cr - Ar 3 )] 



(14) 



with 2c=(£- 1)(£ + 2) and i > 2. 

For perturbations with I > 0, we can show explicitly 
that all the effective potentials V(x) = V{r(x)) are pos- 
itive definite. For scalar perturbations with £ — 0, how- 
ever, the effective potential has one zero point xq and it 
is negative for x > Xq. 



III. NUMERICAL AND SEMI-ANALYTICAL 
APPROACHES 

A. Characteristic integration 

In the work 23] a very simple but at the same 
time very efficient way of dealing with two-dimensional 
d'Alembertians has been set up. Along the general lines 
of the pioneering work [24| , the author introduced light- 
cone variables u = t — x and v = t + x, in terms of 
which all the wave equations introduced have the same 
form. We call V the generic effective potential and ipg the 
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generic field, and the equations can be written, in terms 
of the null coordinates, as 



d 2 



dudv 



ipe(u,v) = V(r(u,v))ipi(u,v) 



(15) 



In the characteristic initial value problem, initial data 
is specified on the two null surfaces u — uo and v = Vq. 
Since the basic aspects of the field decay are independent 
of the initial conditions (this fact is confirmed by our 
simulations), we begin with a Gaussian pulse 01111 = no 
and set the field to zero on v = vq, 



i>i{u = Uo,v) = cxp 



ipi(u, v = v ) = exp 



(.V - Vcf 

2a 2 



2 - 



2a 2 



(16) 



(17) 



Since we do not have analytic solutions to the time- 
dependent wave equation with the effective potentials 
introduced, one approach is to discretize the equation 
(|15fl . and then implement a finite differencing scheme to 
solve it numerically. One possible discretization, used for 
example in 



lcrically. I 



ipe(N) 

- A 2 V(S) 



MW) + ME) 



+ 0(A 4 ), (18) 



while the second time derivative up to 2 nd order is 

tpi(t + At, x) - 2ip t (t, x) + i(ji(t - At, x) 



ipe{t,x) 



At 2 



(20) 

Given ipt(t = to,x) and ipi — (~t = to — At,x) (or 
tpi(t = t ,x)), we can use (IT§|) and ifflty discretization 
to solve (|l(J|l and calculate ipi(t — to + At, x). This is the 
basic algorithm. At each interaction, one can control the 
error by using the invariant integral (the wave energy) 
associate to itTT 



E 



($(f, x) f + (jji(t, x)) 2 + V(x)Mt, x? 



dx . 
(21) 

We make exhaustive analysis of the asymptotic behav- 
ior of the solutions of IjlOjl with initial conditions of the 
form 



exp 



{x - x f 
2a 2 



ipe(0, x) = C exp 



(x - Xi) 



21 



2a{ 



(22) 



(23) 



The results do not depend on the details of initial condi- 
tions. They are compatible with the ones obtained by the 
usual characteristic integration, with the only, and signi- 
ficative, exception of the I = scalar mode. As we will 
see, its asymptotic value depends strongly on the initial 
velocities ipe(0, x). 



where we have used the definitions for the points: N = 
(u + A, v + A), W = (u + A, v), E = (u,v + A) and 
S = {u,v). With the use of expression l(T%|l. the basic 
algorithm will cover the region of interest in the u — v 
plane, using the value of the field at three points in order 
to calculate it at a forth one. 

After the integration is completed, the values 
ipe(u max ,v) and ipi(u,v max ) are extracted, where 
(vmax) is the maximum value of u (v) on the numerical 
grid. Taking sufficiently large u max and v max , we have 
good approximations for the wave function at the event 
and cosmological horizons. 



C. WKB analysis 

Considering the Laplace transformation of the equa- 
tion (fTtH) . one gets the ordinary differential equation 



d 2 ipj(x) 
dx 2 



[s 2 + V{x)] Mx) = 



(24) 



One finds that there is a discrete set of possible values 
to s such that the function ipg, the Laplace transformed 
field, satisfies both boundary conditions, 



lim ip% e s 



1 



(25) 



B. Non-characteristic integration 

It is not difficult to set up a numeric algorithm to solve 
equation l|l(J|l with Cauchy data specified on a t constant 
surface. We used 4 th order in x and 2 nd in t scheme (see, 
for instance, [2{| for an application of this algorithm to 
seismic analysis). The second spatial derivative at a point 
(t, x), up to 4 th order, is given by 

iflfax) = -^^\iP t {t,x + 2Ax)-\Q^t{t,x + Ax) 

+ 30ip e (t, x) - 16^(t, x - Ax) 

+ip e (t, x-2Ax)] , (19) 



lim ^ e- sx = 1 . (26) 

x — > + oo 

By making the formal replacement s = iui, we have the 
usual quasinormal mode boundary conditions. The fre- 
quencies tu (or s) are called quasinormal frequencies. 

The semi-analytic approach used in this work |l5l 
is a very efficient algorithm to calculate the quasinormal 
frequenc ies, which have been applied in a variety of sit- 
uations [l7f. With this method, the quasinormal modes 
are given by 

luI = (V + P) - i [n + ^ (-2F (2) ) V2 (1 + Q) (27) 
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where the quantities P and Q are determined using 



P= - 



V 



(4) 



V 



(2) 



J 



- + a 

4 



1 

288 



V 



(3) 



(2) 



7 + 60a J 



(28) 



-2F n (2) 6912 



2304 
1 

~ 288 



V, 



(4) 



Vr 



(2) 



J 

(6) 



(3) 



1 4 



v; 



(2) 



77+ 188a 2 



-i 2 



67 + 68a 2 



4a 2 



1 

288 



1 

384 

(3) T/ (5) 



K (3)2 K (4) 



v, 



(2)2 



(2) 



19 + 28a 2 



51 + 100a 2 



(29) 



In the equations (|2"?|) to a = n + 1/2 and the 

superscript (i) denote differentiation with respect to x 
of the potential V{x). The potential and its derivatives 
are then calculated in the point xq, where V(x) is an 
extremum. The integer n labels the modes 




-3,.. 



Re(w) > 
Re(w) < 



(30) 



IV. NEAR EXTREME LIMIT 

To characterize the near extreme limit of the 
Schwarzschild-de Sitter geometry, it is convenient to de- 
fine the dimensionless parameter 5: 



5 = \J\ - 9to 2 A . 



(31) 



The limit < 8 <C 1 is the near extreme limit, where 
the horizons are distinct, but very close. In this regime, 
analytical expressions for the frequencies have being cal- 
culated |26T . l27| . For the scalar and electromagnetic fields, 
the quasinormal frequencies are 



A 



3to 2 A 2 



£{£ + !) 



i \ n + 



(32) 

For the axial and polar gravitational fields, the frequen- 
cies are given by 



3m 2 A 2 



(£ + 2)(£-l)- 



i n + 



(33) 
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Figure 1: Decay of the scalar and electromagnetic fields with 
I = 1, and of the scalar, axial and polar gravitational fields, 
with £ = 2, with the SdS geometry approaching the near ex- 
treme limit. The parameters for the geometry are m = 1.0 and 
5 = 0.01,0.1,0.3. 



They can be compared with the numerical and semi- 
analytic methods presented in the previous section. 
Direct calculation of the wave functions confirms that, in 
the near extreme limit, the dynamics is simple, with the 
late time decay of the fields being dominated by quasi- 
normal modes. All the types of perturbation tend to 
coincide near the extreme limit. Besides, as we approach 
the extreme limit, the oscillation period increases and the 
exponential decay rate decreases. These conclusions, il- 
lustrated in figure [T] for t = 1,2, are consistent with the 
presented in |2rll27|. 

By using a non-linear fitting based in a y 2 analysis, it 
is possible to estimate the real and imaginary parts of 
the n = quasinormal mode. These results can be com- 
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pared with the analytical expressions in the near extreme 
cases. In the figure |2 we analyze the dependence of the 
frequencies with I. The accordance is extremely good. 





Figure 3: Asymptotic solutions &(x, i) obtained by non- 
characteristic numerical integration with I = 0, A = 1CP 4 
and m — 1. The curves correspond to different values of C in 
the initial condition 1236 . 



Figure 2: Curve of Re(ujo) 2 with £(£ + 1), for the scalar, elec- 
tromagnetic and gravitational fields, in the near extreme limit. 
The dots are the numerical results and the solid lines repre- 
sent the analytical results. The parameters for the geometry 
are m = 1.0 and 8 — 0.01. 



V. INTERMEDIARY REGION IN PARAMETER 
SPACE 

A. Scalar Field with £ = 

Only scalar perturbations can have zero total angular 
momentum. Solutions of (|1(J|) with I = leads to a con- 
stant tail, as already shown in 13]. This is confirmed 
in figure [21 The novelty here is the dependence of the 
asymptotic value on the tpe(0,x) initial condition. Fig- 
ure 03 reveals the appearance of the constant value <po for 
large t, and its dependence on ipi(Q,x). Note that <fio 
falls below 10~ 7 for ipe(0,x) — 0. These results are in 
accordance with the analytical predictions of 0] > which 
give 

r f Tc ■ ds 
^(oo,r) = ^/ i>(0,s)s—- . (34) 
r c Jo n \s) 



B. Fields with £ > 

We can have scalar and vector fields with angular mo- 
mentum i = 1, and with I > 1, it is possible to introduce 
also gravitational fields. Their behavior is described in 
general by three phases. The first corresponds to the 
quasinormal modes generated from the presence of the 
black hole itself. A little later there is a region of power- 
law decay, which continues indefinitely in an asymptotic 
flat space. In the presence of a positive cosmological con- 
stant, however, an exponential decay takes over in the 
latest period. 



1=1 1 = 2 




K K 

Figure 4: Graphs of the real and imaginary part of the fun- 
damental frequencies (n — 0) with k+. The dotted lines are 
the near extreme results, the dots are the numerical results 
and the continuous curves are the semi-analytic results. In 
the graphs, m = 1.0. 



Some qualitatively different effects show up when we 
turn away from the near extreme limit. For a small cos- 
mological constant the asymptotic behavior is dominated 
by an exponential decaying mode rather than by a quasi- 
normal mode. Such exponential decays will characterize 
the de Sitter space in general. We will later comment 
about the pure de Sitter limit, which is exactly solvable, 
though we do not have definite answers about the respec- 
tive behavior. The exponential decay also characterizes 
the electromagnetic as well as the gravitational pertur- 
bations. 

As the separation of the horizons increases, the quasi- 
normal frequencies deviate from the predicted by expres- 
sions (|32ll and l|33|l . In figure 01 this is illustrated for 
I = 1,2. 

It is interesting to compare the values obtained for 
the fundamental modes using the numerical and semi- 
analytic methods. We find that the agreement between 
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Numerical 


Semi-analytical 


i 


A 


Re(w) 


-Im(oj) 


Re(cj) 


-Im(tj) 


1 


1.000E-5 


2.930E-01 


9.753E-01 


2.911E-01 


9.780E-02 




1.000E-4 


2.928E-01 


9.764E-02 


2.910E-01 


9.797E-02 




1.000E-3 


2.914E-01 


9.726E-02 


2.896E-01 


9.771E-02 




1.000E-2 


2.770E-01 


9.455E-02 


2.753E-01 


9.490E-02 




l.OOOE-1 


8.159E-02 


3.123E-02 


8.144E-02 


3.137E-02 


2 


1.000E-5 


4.840E-01 


9.653E-02 


4.832E-01 


9.680E-02 




1.000E-4 


4.833E-01 


8.948E-02 


4.830E-01 


9.677E-02 




1.000E-3 


4.816E-01 


8.998E-02 


4.809E-01 


9.643E-02 




1.000E-2 


4.598E-01 


8.880E-02 


4.592E-01 


9.290E-02 




l.OOOE-1 


1.466E-01 


3.068E-02 


1.466E-01 


3.070E-02 


3 


1.000E-5 


6.769E-01 


8.662E-02 


6.752E-01 


9.651E-02 




1.000E-4 


6.754E-01 


8.654E-02 


6.749E-01 


9.647E-02 




1.000E-3 


6.732E-01 


8.660E-02 


6.720E-01 


9.611E-02 




1.000E-2 


6.437E-02 


9.200E-02 


6.428E-02 


9.235E-02 




l.OOOE-1 


2.091E-02 


3.054E-02 


2.091E-02 


3.056E-02 



Table I: Fundamental frequency (n = 0) of the scalar field, 
obtained using the numerical and semi-analytical methods. In 
this table, m = 1.0. 





Numerical 


Semi-analytical 


i 


A 


Re(w) 


-Im(cj) 


Re(uj) 


-Im(tj) 


1 


1.000E-5 


2.481E-01 


9.226E-02 


2.459E-01 


9.310E-02 




1.000E-4 


2.481E-01 


9.223E-02 


2.457E-01 


9.307E-02 




1.000E-3 


2.475E-01 


9.176E-02 


2.448E-01 


9.270E-02 




1.000E-2 


2.374E-01 


8.839E-02 


2.352E-01 


8.896E-02 




l.OOOE-1 


8.035E-02 


3.027E-02 


8.023E-02 


3.033E-02 


2 


1.000E-5 


4.577E-01 


8.985E-02 


4.571E-01 


9.506E-02 




1.000E-4 


4.575E-01 


8.991E-02 


4.569E-01 


9.502E-02 




1.000E-3 


4.559E-01 


9.439E-02 


4.551E-01 


9.464E-02 




1.000E-2 


4.371E-01 


8.941E-02 


4.364E-01 


9.074E-02 




l.OOOE-1 


1.458E-01 


3.037E-02 


1.458E-01 


3.038E-02 


3 


1.000E-5 


6.578E-01 


8.365E-02 


6.567E-01 


9.563E-02 




1.000E-4 


6.576E-01 


8.349E-02 


6.564E-01 


9.559E-02 




1.000E-3 


6.547E-01 


8.399E-02 


6.538E-01 


9.520E-02 




1.000E-2 


6.276E-02 


8.852E-02 


6.267E-01 


9.125E-02 




l.OOOE-1 


2.085E-02 


3.039E-03 


2.085E-01 


3.040E-02 



Table II: Fundamental frequency (n — 0) of the electromag- 
netic field, obtained using the numerical and semi-analytical 
methods. In this table, m — 1.0. 



then is good, for the hole range of A. The difference is 
better for the first values of i. This is expected, since the 
numerical calculations are better in this region. In the 
tabled we illustrate these observations with a few values 
of A. It is important to mention that quasinormal fre- 
quencies for the SdS black hole were already calculated 
in a recent paper llfll, applying a variation of the WKB 
method used here |l8j . There are earlier papers calcu- 
lating quasinormal modes in this geometry, for example 

The first higher n modes cannot be obtained from the 
numerical solution, but can be calculated by the semi- 
analytical method. As the cosmological constant de- 
creases, the real and imaginary parts of the frequencies 
increase, up to the limit where the geometry is asymp- 
totically flat. The behavior of the modes is illustrated 
in figure |SJ The behavior of the electromagnetic field is 
similar. 

A x 2 analysis of the data presented in figure shows 
that the massless scalar, electromagnetic and gravita- 
tional perturbations in SdS geometry behave as 





Numerical 


Semi-analytical 


t 


A 


Re(w) 


-Im(w) 


Re(u) 


-Im(a;) 


2 


1.000E-5 


3.738E-01 


8.883E-02 


3.731E-01 


8.921E-02 




1.000E-4 


3.737E-01 


8.880E-02 


3.730E-01 


8.918E-02 




1.000E-3 


3.721E-01 


8.850E-02 


3.715E-01 


8.888E-02 




1.000E-2 


3.566E-01 


8.538E-02 


3.560E-01 


8.572E-02 




1.000E-1 


1.179E-01 


3.020E-02 


1.179E-01 


3.023E-02 


3 


1.000E-5 


5.999E-01 


8.677E-02 


5.992E-01 


9.272E-02 




1.000E-4 


5.996E-01 


8.676E-02 


5.990E-01 


9.269E-02 




1.000E-3 


5.972E-01 


8.971E-02 


5.966E-01 


9.234E-02 




1.000E-2 


5.725E-01 


8.695E-02 


5.718E-01 


8.874E-02 




1.000E-1 


1.900E-01 


3.030E-02 


1.900E-02 


3.032E-02 


4 


1.000E-5 


8.106E-01 


8.810E-02 


8.091E-01 


9.417E-02 




1.000E-4 


8.102E-01 


8.781E-02 


8.087E-01 


9.413E-02 




1.000E-3 


8.070E-01 


8.799E-02 


8.055E-01 


9.376E-02 




1.000E-2 


7.733E-02 


8.714E-02 


7.720E-01 


9.000E-02 




1.000E-1 


2.564E-02 


3.034E-03 


2.563E-01 


3.036E-02 



with t 



(35) 



Table III: Fundamental frequency (n = 0) of the axial and 
polar gravitational fields, obtained using the numerical and 
semi-analytical methods. In this table, m = 1.0. 



,po 



with t — > oo 



with t — > oo 



e fe ^p* with t 



(36) 
(37) 
(38) 



for t sufficiently large. At the event and the cosmological 
horizons t is substituted, respectively by v and u. 



The numerical simulations developed in the present 
work reveal an interesting transition between oscillatory 
modes and exponentially decaying modes. As shown in 
figure [7| as the cosmological constant decreases, the ab- 
solute value of — Im(w) decreases. 

Above a certain critical value of A we do not observe 
the exponential tail, since the coefficient k exp is larger 
than — Im(o;o) thus the decaying quasinormal mode domi- 
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Figure 5: Quasinormal modes of the scalar, axial and polar 
gravitational fields, for higher modes. The parameter for the 
curves are m = 1 and A = 10 -3 . 
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Figure 6: Exponential tail for the scalar and electromagnetic 
fields with 1=1, and for the axial and polar gravitational field 
with i = 2. In the graphs, m = 1.0. 




0.009 0.02 o 

A 



0.002 0.004 0.0005 0.001 

A A 



Figure 7: Approaching of the constants —Im(ujo) and k £xp , 
for the scalar, electromagnetic and gravitational fields, in the 
SdS geometry. Above a certain critical value of A (roughly 
1.7 x 10~ 2 , 4.0 x 10" 3 and 1.2 x 10" 3 respectively, for the 
parameters taken in graphs), a tail it is not observed. For all 
curves, the mass is set to m = 1.0. 
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nates. But for A smaller than this critical value, — Im(wo) 
turns out to be larger than k exp and the exponential tail 
dominates. Certainly, for a small enough cosmological 
constant the exponential tail dominates in the various 
cases considered here. 

Another aspect worth mentioning in the intermediate 
region is the dependence of the parameters k^ p , kf xp , 
fcg^p and fcf° p with £ and k c . The results suggest that 
the k exp are at least second diffcrcntiable functions of k c . 
Therefore, close to k c = 0, we approximate 

K c xp (k c )^1(k c + c sc k 2 c ) , (39) 



Figure 8: Dependence of k exp with k c and i, in the SdS ge- 
ometry. The symbols indicate the numerical values, and the 
solid lines are the appropriate fittings. For the left graph: 
kll p = 1.077 x 10~ 4 + 0.984«: c + 3.545k 2 , kf xv = -1.863 x 
10" 4 +2.010k c + 2.608«: 2 andkf xp = -1.959x 10 _4 +3.028fi: c + 
2.978k 2 . For the center graph: kf xp = 2.082 x 10" 4 + 
1.988k c + 6.141k 2 , kf xp = 2.712 x 10 -4 + 2.974k c + 8.284k 2 
and kf xp = 3.737 x 10" 4 + 3.984k c + 4.517k 2 . For the 
right graph: k^ p = 2.616 x 10" 4 + 2.974k c + 9.895k 2 and 
Kxp = 4.484 x 10" 4 + 3.896k c + 18.92k 2 . In the graphs, 
m = 1.0. 



(40) Previous results are illustrated in figure |SJ 
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Figure 9: Approaching the asymptotically flat geometry. 
Straight lines in the log-log graphs indicate power-law decay. 
In the graphs, m = 1.0. 



VI. APPROACHING THE ASYMPTOTICALLY 
FLAT GEOMETRY 

Scalar fields in the SdS geometry near the asymptoti- 
cally flat limit wore studied in [lj, 13 • In this case there 
is a clear separation between the event and the cosmo- 
logical horizons, such that 

6 = r J—I± > 50 . (41) 

r+ 

A new qualitative change occurs in this regime, namely 
a decaying phase with a power-law behavior. Such a 
phase occurs between the quasinormal mode decay and 
the exponential decay phases. The field cannot be simply 
described by a superposition of the various modes, which 
would imply a domination of the power-law phase. This 
is illustrated in figure |5J 



VII. CONCLUSIONS 

We have identified three regimes, according to the 
value of A for the decay of the scalar, electromagnetic 
and gravitational perturbations. Near the extreme limit 
(high A), we have analytic expressions for the effective 
potentials and to the quasinormal frequencies. The de- 
cay is entirely dominated by the quasinormal modes, that 
is, oscillatory decaying characterized by a non vanishing 
real part of the quasinormal frequency. 

In an intermediary parameter region (lower A), the 
wave functions have an important qualitative change, 
with the appearance of an exponential tail. This tail 
dominates the decay for large time. Near the asymptot- 
ically flat limit (A <C 1), we see an intermediary phase 
between the quasinormal modes and the exponential tail 
— a region of power-law decay. When A = 0, this region 
entirely dominates the late time behavior. 

Finally, for scalar fields with I = a constant decay 
mode appears, and its value depends on the tpe(0,x) ini- 
tial condition. Figure[21reveals the appearance of the con- 
stant value 4>q for large t, and its dependence on ifit(0, x). 
The value of 0o falls below 1CU 7 for tpe(0,x) = 0. 

An immediate extension of the work developed in this 
paper is the in treatment of field propagation in the exte- 
rior of charged and asymptotically de Sitter black holes. 
In this case, the electromagnetic and gravitational per- 
turbations are necessarily coupled. It would be interest- 
ing to see if the general picture presented here is still 
valid in this more general context. 

Acknowledgments 

This work was supported by Fundacao de Amparo 
a Pesquisa do Estado de Sao Paulo (FAPESP) and by 
Conselho Nacional de Desenvolvimento Cientffico e Tec- 
nologico (CNPq), Brazil. 



[1] T. Regge and J. A. Wheeler, Phys. Rev. 108, 1063 (1957). 

[2] S. Chandrasekhar, The Mathematical Theory of Black 
Holes, (Oxford University Press, Oxforf, 1983). 

[3] K. D. Kokkotas and B. G. Schmidt, Living Rev. Relativ- 
ity 2 (1999), gr-qc/9909058 

[4] A. Riess et al. Astrophy. J. 116, 1009 (1998), 
astro-ph/9805201 S. Perlmutter et al., Astrophys. J. 
517, 565 (1999),| astro-ph/9812133l 

[5] V. Car doso and J. P. S . Lemos, Phys. Rev. D 64, 084017 
(2001), gr-qc/0105103 

[6] B. Wang, (J. Lin and E. A bdalla, Phys. Lett. B, 481, 
79 (2000), hep-th/0003295 R. A. Konop lya, Phys. Rev. 
D 66, 044009 (2002), hep-th/0205142 E. Berti and 
K. D. Kokkot as, Phys. Rev. D 67, 064020 (2003), 
|gr-qc/0301052l V. Cardoso, R. Konopl ya and J. P. S. 
Lemos, Phys. Rev. D 68, 044024 (2003), gr-qc/0305037 

[7] V. Cardoso and J. P. S. Lemos, Class. Quant. Grav. 



18, 5257 (2001), gr-qc/0107098 V. Card oso and J. P. S 
Lemos, Phys. Rev. D 63, 124015 (2001), gr-qc/0101052 
J. M. Zhu, B. Wang and E. Abdalla, Phys. Rev. D 
63, 124004 (2001), |hep-th/010iT33l B. Wang, E. Ab- 
dalla and R. B. Mann, Phys. Rev. D 65, 084006 (2002), 
|hep-th/0107243| 
[8] B. Wang, C. Molina and E. Ab dalla, Phys. Rev. D 63, 

084001 (2001), hep-th/0005143 
[9] H. Nariai, Sci. Rep. Tohoku Univ. 34, 160 (1950); H. 
Nariai, Sci. Rep. Tohoku Univ. 35, 62 (1951). 

[10] E. W. Leaver, Phys. Rev. D 34, 384 (1986). 

[11] E. S. C. Ching, P. T. Leung, W. M. Suen and K. Young, 
Phys. Rev. D 52, 2118 (1995), gr-qc/9507035 

[12] G. Gamow, Zeits. F. Phys. 51 204 (1928). 

[13] P. R. Brady, C. M. Chambers, W. Krivan and P. Laguna, 
Phys. Rev. D 55, 7538 (1997), gr-qc/9611056 

[14] P. R. Brady, C. M. Chambers, W. G. Laarakkers and E. 



9 



Poisson, Phys. Rev. D 60, 064003 (1999), gr-qc/9902010 
[15] B. F. Schutz and C. M. Will, Astrophys. J. 291, L33 
(1985). 

[16] S. Iyer and C. M. Will, Phys. Rev. D 35, 3621 (1987). 
[17] S. Iyer, Phys. Rev. D 35, 3632 (1987); K. D. Kokkotas 

and B. F. Schutz, Phys. Rev. D 37, 3378 (1988). 
[181 R. A. Konopl ya, Phys. Rev. D 68, 024018 (2003), 

gr-qc/0303052 

[19] A. Zhidenko, Quasi-normal modes of Schwarzschild-de 
Sitter black holes, gr-qc/0307012 

[20] S. Yoshida and T. Futamase, Numerical analysis of quasi- 
normal modes in nearly extremal Schwarzschild-de Sitter 
spacetimes, gr-qc/0308077 

[21] E. Abdalla, B. Wang, A. Lima-Santos and W.G. Qiu, 
Phys. Lett. B 538, 435 (2002), hep-th/0204030 E. Ab- 
dalla, K.H.C. Castello-Branco, A. Lima-Santos, Mod. 



Phys. Lett. A 18, 1435 (2003), gr-qc/0301130 
[22] R. Ruffini, J. Tiomno and C. Vishveshwara, Lettere al 

Nuovo Cimento 3, 211 (1972). 
[23] C. Gundlach, R. Price and J. Pullin, Phys. Rev. D 49, 

883 (1994). 
[24] R. Price, Phys. Rev. D 5, 2419 (1972). 
[25] A.R. Levander, Geophysics 53, 1425 (1988). 
[26] V. Car doso and J. P. S . Lemos, Phys. Rev. D 67, 084020 

(2003), gr-qc/0301078 
[271 C. Molina, Phys Rev. D 68, 064007 (2003), 

gr-qc/0304053 

[28] H. Otsuki and T. Futamase. Prog. Theor. Phys. 85. 771 
(1991); I. Moss and J. Norman, Class. Quant. Grav. 19, 
2323 (2002). 



